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1 . 0 INTRODUCTION 


1.1  Technical  Background 

For  a comprehensive  study  of  the  ionized  state  of 
the  upper  atmosphere,  it  is  necessary  to  consider  all  rele- 
vant fundamental  processes  which  are  responsible  for  main- 
taining that  state  of  ionization.  As  the  earth’s  magnetic 
field  serves  to  confine  the  plasma,  the  electron  density  and 
temperature  profiles  exhibit  marked  differences  for  different 
magnetic  latitudes.  In  consequence,  to  understand  the  dynam- 
ic processes  which  involve  transport  of  electrons  from  one 
region  of  the  earth  to  another  region,  one  must  make  a dis- 
tinction among  three  regions  of  the  atmosphere.  These  are: 

a)  Equatorward  of  the  precipitation  zone. 

b)  The  precipitation  region. 

c)  The  polar  region. 

The  distinguishing  features  of  these  regions  are  related  to 
the  configurations  of  the  magnetic  fields.  For  the  equator- 
ward  side  of  the  precipitation  zone,  the  field  lines  below 
the  plasma  pause  are  essentially  closed.  In  contrast,  the 
field  lines  in  the  precipitation  zone  are  presumably  con- 
nected to  the  plasma  sheet,  while,  in  the  polar  region  the 
field  lines  may  find  a possible  connection  either  with  the 
magnetosphere  or  the  interplanetary  magnetic  fields.  One  of 
the  careful  considerations  engendered  by  these  differing  con- 
figurations in  dealing  with  transport  of  electrons  is  the 
corresponding  differing  boundary  conditions  and  the  associated 
problems  of  the  energy  balance.  A study  and  analysis  of  the 
transi>ort  phenomenon  and  the  morphology  of  the  upper  atmo- 
sphere was  carried  out  in  cooperation  with  the  Electrical 
Processes  Branch  of  AFGL.  In  the  following  few  sections  we 


1 


give  a brief  exposition  of  some  of  the  basic  ideas  that  are 
involved  in  this  kind  of  study  and  also  sketch  our  method  of 
approach. 


1.2  The  Basic  Problem 

For  a given  steady  state  flux  of  hyperthermal  elec- 
trons whose  initial  energy  spectrum  is  prescribed,  if  we  as- 
sume that  the  neutral  species  and  ambient  plasma  profiles  are 
also  given  (or  prescribed),  then  what  is  the  final  steady 
state  of  the  hyperthermal  electrons  as  they  travel  along  the 
magnetic  field  lines?  We  allow  the  hyperthermal  electrons  to 
interact  with  the  neutral  species  and  the  ambient  plasma,  but 
otherwise  they  are  correlationless,  i.e.  hyperthermal  elec- 
trons do  not  interact  with  each  other. 

[ 1.3  The  Basic  Quantity  Sought 

The  basic  quantity  we  wish  to  find  is  the  electron 
distribution  function,  which  is,  the  expected  number  of  elec- 
trons per  cm*  per  electron  volt  per  steradian.  Once  the  dis- 
tribution function  is  determined,  then  any  other  macroscopic 
quantity  is  determined  by  averaging  with  respect  to  that. 

1.4  Some  Problems  of  Interest  Addressed 

There  are  several  relevant  problems  associated  with 
the  three  regions  of  the  upper  atmosphere  mentioned  before. 

The  basic  question  pertains  to  the  energy  balance  of  the  spec- 
tral (energy)  flux  of  the  hyperthermal  electrons  associated 
with  the  energy  budget  of  the  neutral  and  the  thermal  cunbient 
plasma  environment.  By  that  we  mean  that  any  eunount  of  energy 
imparted  to  the  environment  by  the  hyperthermal  electrons  must 
be  accountable.  This  energy  may  appear  either  in  the  form  of 
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ambient  plasma  heating  or  excitation  (also  possible  ioniza- 
tion) of  neutral  species  with  a possible  subsequent  optical 
emission.  To  elaborate  on  this  point  let  us  first  consider 
the  role  played  by  the  hyperthermal  electrons  produced  by 

photoionization.  About  one-third  of  the  solar  energy  spec- 

o o 

trum  between  soft  x-rays  (~14A)  and  EUV  (-1100A)  is  expended 


into  the  photoionization  and  (or)  excitation  of  the  atmo- 
spheric neutral  species.  In  general,  the  energy  range  of  the 
photoelectron  fluxes  lies  between  zero  and  200  eV  (electron 
volts)  and  the  significant  altitude  range  where  they  are  pro- 
duced is  between  about  120  km  to  about  400  km.  In  a real 
atmosphere,  as  these  photo electrons  travel  along  the  magnetic 
field  lines,  they  lose  their  energy  by  collisions  with  neutral 
species  and  also  by  Coulomb  interaction  with  the  relatively 
cooler  ambient  plasma.  The  collisions  with  neutral  species 
may  result  in  excitation  (and)  or  further  ionization.  If  it 
is  an  excitational  collision  then  the  photoelectrons  will 
lose  their  energy  in  discrete  jumps  determined  by  the  excita- 
tion threshold  energy.  The  excitation-energy  may  subsequently 

O 

be  lost  by  either  optical  emission  (such  as  6300A  red  line  of 
state)  or  be  quenched  by  collisional  deexcitation  (such  as 
Nj  vibrational  states).  If  it  is  an  ionizational  collision 
then,  the  total  energy  of  emerging  electrons  (which  may  have 
continuous  energy  distribution)  plus  whatever  goes  into  ex- 
citation-energy  of  the  ion,  must  balance  the  energy  of  the 
impacting  electron.  In  consequence,  there  are  several  chan- 
nels into  which  an  impacting  electron  can  distribute  its  en- 
ergy in  interaction  with  neutral  species.  With  respect  to 
the  ene-’gy  loss  to  the  cooler  ambient  plasma,  the  loss  process 
is  continuous.  However,  the  important  point  is  that  the 
Coulomb  interaction  will  result  in  raising  the  ambient  elec- 
tron temperature  more  efficiently  than  that  of  the  ions.  As 
an  specific  example  let  us  consider  the  transport  of  photo- 
electro is  from  the  sunlit  side  to  the  conjugate  dark  side. 
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I Some  of  the  data  analyzed  from  the  electrostatic  analyzer  on- 

( board  the  INJUNS  satellite  exhibited  that  the  hyperthermal 

electron  fluxes,  in  the  altitude  range  of  1000  Jon  to  2500  km, 
with  energies  greater  than  2 eV  were  a consequence  of  the 
photoelectrons  which  were  transported  from  the  sunlit  side  of 
the  earth.  Thus,  in  the  accessed,  dark-regions  of  the  top- 
side ionosphere,  the  observed  temperatures  were  between  3000 
to  4000® K with  highest  temperatures  near  the  trough.  The  un- 
accessed regions  were  essentially  isothermal  with  the  tempera- 
ture in  the  neighborhood  of  2000®K.  As  a coauthor  to  W.  J. 
Burke  and  R.  C.  Sagalyn,  this  study  of  observations  has  been 
submitted  as  a paper  to  J.  Geophys.  Res.  (1977).  Another 
interesting  exaunple  of  the  transport  phenomenon  is  the  excita- 
tion of  the  stable  auroral  arcs  (SAR).  As  is  well  known, 
during  high  Kp-indices  the  ring  current  protons  move  inward 
into  the  cool  plasmasphere  and  the  ring  current  is  eroded  by 
the  cyclotron  wave  turbulence  instability.  The  resulting 
transfer  of  energy  to  plasmaspheric  electrons  is  conducted 
down  into  the  lower  ionosphere  where  the  ambient  electrons 
are  heated  to  as  much  as  5000°K.  If  the  energy  density  of 
the  ambient  electrons  is  sufficiently  high  (about  0.02  ergs) 
the  high  energy  tail  of  the  distribution  will  then  contain 

sufficient  number  of  electrons  (about  2 eV)  to  excite  the  *0_ 
o ^ 

(6300A)  line  of  the  atomic  oxygen.  Of  course,  the  reaction 

rate  will  directly  depend  on  the  density  of  the  atomic  oxygen. 
For  that  reason  SAR  arcs  are  observed  in  the  altitude  range 
of  about  250-600  km.  However,  the  more  important  point  is 
that  here  the  transport  phenomenon  combines  with  the  heat  con- 
duction and  the  energy  source  is  the  ring  current.  A paper 
on  this  subject  has  been  submitted  to  J.  Geophys.  Res.  (1977). 
To  mention  two  other  examples  of  interest,  one  is  the  second- 
ary production  of  electrons  and  their  distribution  in  the  pre- 
cipitation zone  and  the  other  is  the  escape  of  plasma  in  the 
polar  region.  In  the  precipitation  zone,  if  we  assume  that 
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the  flux  of  high  energy  electrons  (1-50  keV)  is  of  known  na- 
ture, then  this  boundary  condition  is  sufficient  to  determine 
the  distribution  (and  secondary  electron  production  rate) 
along  the  entire  path.  Here  one  can  envisage  an  interesting 
possibility  of  an  inverse  problem,  i.e.  given  the  results  of 
some  rocket  experiment  in  the  auroral  (precipitation)  zone, 
can  one  determine  the  nature  of  the  source  of  electrons?  We 
investigated  the  possibility  of  such  a type  of  data  profile 
inversion. 


2.0  ELECTRON  TRANSPORT  EQUATION 


We  now  give  a brief  sketch  of  the  derivation  of  the 
transport  equation.  From  any  closed  volume  fixed  in  space 
the  electrons  flow  out,  whereby  the  advective  term  equal  to 
the  divergence  of  the  current  contributes  to  the  loss  of  elec- 
trons in  the  phase  space. 

The  electron  transport  equation  emerges  when  we 
balance  the  source  and  sink  terms.  Thus,  we  have  the  contin- 
uity equation  which  is  given  by 

Div.  of  Electron  Current  + Electron  Sinks  = 

Source  Function  + SES  + PES  (1) 

where  SES  represents  the  secondary  electron  spectrum  and  PES 
the  primary  electron  spectrum.  We  now  define  the  photoelec- 
tron distribution  (PED)  as 

PED  = expected  number  of  electrons  (excluding 

ambient  plasma  electrons)  per  unit  volume 
per  electron  volt  per  steradian.  (2) 

Let  f(r,  E,  y)  denote  PED,  where  r is  the  position  vector,  E 
the  energy  and  y the  pitch  angle  of  the  electron.  Various 
components  of  Eq.  (1)  are  then,  related  to  f(r,  E,  y)  as 
follows : 


Div.  of  Electron  Current  = V«ufCr,  E,  y)  (3) 

where  u is  the  velocity  vector  of  the  electron  and  V is  the 
spatial  gradient  operator.  Let  denote  the  total  collision 
rate  of  electrons  with  the  medium.  Then 


Electron  Sinks  = v^f(r,  E,  y) 


(4) 


If  nj  represents  the  number  density  (no/cm®)  of  neutral 

(ex)  (ion) 

species,  Qjj(E)  the  cross-section  for  excitation  and  Qjj(E) 

the  ionization  cross-section  corresponding  to  J-^  state  of 


(e) 

species  I,  Qj(E)  the  elastic  scattering  cross-section  of 
species  I and  r(E)  (=  AE/At)  the  rate  of  energy  loss  to  the 
ambient  plasma,  then  the  total  collision  rate  is  given  by 


Let  R(E',  fl'  E,  n)  generally  represent  the  proba- 
bility (or  the  redistribution  function)  that  an  electron  of 
initial  energy  E'  and  momentum  unit  vector  upon  collision 
emerges  with  energy  E and  momentum  unit  vector  fi.  Then  the 
source  function  (emission  term)  in  Eq.  (1)  is  given  by 

00 

SOURCE  FUNCTION  = / dE'  / dfl'  R(E’,  O'  -*■  E,  «) 

HTT  Q 


f(r,  E' , 0)  (6) 

where  m (=  9.1  x lo”^®  gm)  is  the  mass  of  the  electron. 

For  elastic  scattering  from  neutral  species 
R(E',  ft'  E,  ft)  = 6(E'-E)  Rg(fi'  «)  (7) 

where  ^)  is  the  angular  phase  function  corresponding 

to  the  elastic  scattering. 

For  scattering  by  plasma  we  assume  that  energy  lost 
by  the  incident  electron  is  AE.  Then 

R(E',  ft'  E,  ft)  = 6(E'-E-AE)  Rp(fl*  ft)  (8) 

where  Rp(ft'  -►  ft)  is  the  angular  phase  function  corresponding 
to  scattering  by  ambient  plasma. 
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For  an  inelastic  collision  in  which  the  neutral 
species  is  left  in  an  excited  state, 

R(E',  E,  a)  = 6(E*-E-Wjj)  fi)  (9) 

where  Wjj  is  the  threshold  energy  and  R^j  is  the  corresponding 
angular  phase  function. 

Substituting  the  redistribution  functions  for  the 
various  processes  as  given  by  Eqs.  (7)  - C9)  into  Eq.  (6)  we 
obtain 

SOURCE  FUNCTION  = / dO*  R^Cft'  fli)  I n^Qj^fCr,  E,  O') 

, , /STE-T-g— T (ex) 

+ ^ / dfi'  Rp(«'  -*■  £)  f(r,  E + AE,  Q')  . (10) 

To  calculate  the  secondary  electron  spectrum  let  us  write  the 
redistribution  function  as 

R(E',  ft'  -*■  E,  ft)  = a(E'  E)  R^(ft'  -*■  ft)  (11) 

where  o(E'  E)  is  differential  cross-section  (cf.  Opal  et  al, 
1971)  for  producing  a secondary  electron  of  energy  E by  an 
electron  of  initial  energy  E'.  R (ft'  -►  ft)  is  the  corresponding 
angular  phase  function  for  the  given  process.  We  assume  that 
in  the  primary  electron  spectrum  (PES)  the  maximum  energy  of 
an  electron  is  E^^.  If  an  electron  of  initial  energy  E'  has  an 
ionizing  collision  with  a neutral  species  with  two  electrons 
of  energies  E^  and  Eg  exiting,  and  ionized  species  left  in 
some  excited  state,  then  the  conservation  of  energy  states 
that 

E'  = E + E„  + W,T  (12) 

P b XJ 

where  Wjj  is  the  ionization  threshold  plus  the  energy  corre- 
sponding to  the  excited  state.  We  now  wish  to  know  how  many 
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secondary  electrons  of  energy  Ep  are  produced  per  cm*  per  sec 
per  electron  volt  by  electrons  of  energy  E'.  The  value  of  E* 
may  range  from  Ej^j  to  some  lowest  value  permitted  by  the  energy 
conservation  equation  (12).  Evidently  the  second  electron  of 
energy  Eg  may  end  up  in  any  part  of  the  energy  spectrum 
bounded  by 

0 < Eg  < Ej^  - Ep  - Wjj  . (13) 

This  implies  that  in  order  for  an  ionization  collision  to 
occur  E'  must  obey  the  inequality 

Ep  * W„  < E'  < . (1-) 


Consequently,  the  total  number  of  electrons  of  energy  Ep  pro- 
duced per  cm^  per  sec  per  electron  volt  per  steradian  are 


Ew— Ep— Wt  j E^  __ 

M P IJ  M /TrT 

SES  = ^ I n,  / dft*  R„(n'  il)  / dE„  / dE*  a(E 

IJ  ^ 0 ^ 


’-Ep) 


6(E'  - Ep  - Eg  - Wjj)  f(r,  E',  fl') 


where  Dirac's  delta  function  6(E'  - Ep  - Eg  - W^j)  ensures  the 
conservation  of  energy.  Integration  in  (15)  with  respect  to 
E'  is  rather  trivial.  With  the  change  of  variable  Ep  + Eg  = E" 
we  get 

= h J,  "1/““'  L “s'e’  - «)  s/s  <E"  ♦ Wjj 

I ,u  Ep 


• o(E"  + Wjj  -►  Ep)  f(r,  E"  + Wjj,  fl ' ) . (16) 

Now  substituting  the  explicit  expressions  for  vari- 
ous constituents  of  the  continuity  Eq.  (1)  as  given  by  (3), 
(4),  (10)  and  (lo)  we  get  the  following  electron  transport 
equation 
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V • u f(r,  E,  a)  + V.  f(r,  E,  ft)  = i-  / dft*  R (ft*-ft)  \/p  I n^Q^CE) 

— — “ T — — HIT  ' e—  — 

• f(r,  E,  ft') 

-1  py (ex) 

♦ 5?  / I,  "I^IJ  <E*W„)  R„cn'-.a>  f(r.  Etw.j,  B') 

I,  J 

+ / dft'  ^ 

E„-Wjj  

+ nj/dft'  / dE'  RgCn'-fft)  \/|(E'+Wj.j)a(E'+Wjj-*-E)f(r,E'+Wj.j,n’) 

+ P(E,  r)  (17) 

where  P(E,  r)  denotes  PES.  A simpler  version  of  this  equation 
has  been  discussed  in  the  literature  (cf.  McCormick  and  Kuscer 
(1965),  Case  and  Zweifel  (1967),  Mika  (1961),  Zelazny  et  al 
(1961)  and  Kanal  (1970)). 

In  writing  down  Eq.  (17),  we  have  neglected  the 
roles  of  electric  field  (e)  and  magnetic  field  (B).  In  gen- 
erals, however,  one  must  add  the  Lorentz  force  term 

F = q(e  + u B)  (18) 

to  the  left  hand  side  of  Eq.  (17)  so  that 

V • u f(r,  E,  ft)  + V.  f(r,  E,  ft)  + i F*V..  f(r,  E,  ft) 


- R„(r,  E,  ft) 


(19) 


where  is  velocity  gradient  operator  and  R^  represents  the 
entire  right  hand  side  of  Eq.  (17). 


3.0  COMPLETED  TASKS 


For  a pitch  angle  scattering  of  electrons,  which  is 
linearly  anisotropic,  we  have  complete  solutions  of  the  trans- 
port equation  (17)  for  all  three  regions  of  the  upper  atmo- 
sphere, i.e.  the  temperate  zone,  the  precipitation  zone  and 
the  polar  region.  A computer  program  is  available,  which  is 
adapted  to  the  AFGL  CDC  6600  computer.  It  was  our  original 
intention  to  run  this  program  completely  without  any  bugs  and 
make  a permanent  file  for  AFGL's  use.  However,  it  was  felt 
that  for  a four  month's  effort  it  would  be  more  useful  to  con- 
struct complete  theories  first  for  the  precipitation  and  the 
polar  regions.  This  decision  was  dictated  by  the  existence 
of  the  wealth  of  INJUNS  data  and  the  anticipated  data  on  elec- 
tric fields  and  plasma  bulk  motion  from  R.  C.  Sagalyn's  exper- 
iments on  S3-2  and  S3-3  Air  Force  Satellites.  Furthermore,  it 
was  also  felt  that  our  model  of  pitch  angle  scattering  which 
is,  linearly  anisotropic,  was  not  adequate  for  the  entire  use- 
ful energy  spectrum  of  hyperthermal  electrons.  For  instance, 
for  high  energy  electrons  (>50  eV)  the  angular  differential 
cross-sections  are  highly  peaked  in  the  forward  direction. 

This  is  clearly  seen  in  the  experiments  of  Crooks  and  Rudd 
(1972)  in  which  electrons  of  energy  between  50  eV  and  800  eV 
were  bombarded  on  helium  gas.  Their  results  are  shown  in 
Figures  1 to  7 . If  we  consider  the  elastic  scattering  first, 
we  see  that  for  all  energy  ranges  (50  eV  to  200  eV)  the  dif- 
ferent;.al  cross-sections,  as  shown  in  Figure  1,  assume  the 
largest  values  for  scattering  angles  close  to  0° . Here  the 
scattering  angle  0 is  measured  with  respect  to  the  initial 
direction  of  the  electron  velocity  vectors  as  shown  in  Figure 
1.  In  particular,  for  50  eV  electTOns  the  probability  for 
forward  scattering  is  approximately  20  times  larger  than  the 
probability  for  backscatter.  This  factor,  of  course,  increases 
with  increasing  energy  of  the  incident  electron.  For  an  ex- 
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I citation  collision.  Figure  2 shows  the  differential  cross- 

; sections  for  50  eV  electrons  producing  2^S,  2*S,  2^P  and  2*P 

states  of  the  helium  atom.  The  important  fact  to  note  is 

I 

that,  indeed,  the  cross-sections  are  forward  peaked  (8  = 0®), 

I but  assume  a lower  value  in  the  lateral  direction  than  the 

backscatter  direction.  Figure  3 shows  the  similar  tendency 
for  100  eV  electrons.  Figures  4 to  6 illustrate  the  nature 
of  the  angular  dependence  of  doubly  differential  cross-section 
for  an  ionizing  collision.  For  the  seune  species  (He) , we  pre- 
sent the  angular  differential  cross-sections  as  obtained  ex- 
perimentally by  Peterson  et  al  (1971,  1972).  These  are  shown 
in  Figures  8 and  9.  There  exists  a considerable  amount  of 
‘ experimental  data  (c.f.  Lassettre  et  al,  1964a,  b,  1968)  which 

I can  be  used  for  constructing  a realistic  angular  model  of 

' scattering  for  the  purpose  of  calculating  the  photoelectron 

I fluxes  in  the  upper  atmosphere. 

The  assumption  that  the  photo electro ns  below  about 
: 10  eV  are  isotropically  scattered  by  neutrals  is  fairly  accu- 

; rate.  However,  this  assumption  is  not  valid  for  the  Coulomb 

interaction.  When  the  energy  of  an  electron  has  degraded  be- 
low about  15  eV,  the  Coulomb  interactions  with  ambient  thermal 
plasma  begin  to  play  an  important  role.  For  higher  energy 
electrons  (about  about  15  eV),  Spitzer's  (1956)  model  of  small 
emgle  scattering  (by  plasma)  may  be  deemed  satisfactory. 

Along  the  same  line,  further  developments  by  Butler  and 
Buckingham  (1962)  and  the  empirical  model  by  Swartz  et  al 
(1971)  for  the  rate  of  energy  loss  (see  Figure  10)  art  indis- 
pensable for  modelling.  However,  below  10  eV  at  close..'  exam- 
ination of  the  pitch  angle  scattering  by  the  thermal  plasma 
is  essential.  In  view  of  the  fact  that  the  maximum  heating 
rate  of  the  ambient  plasma  by  hyperthermal  electrons  occurs 
in  the  neighborhood  of  2 to  3 eV  (except  for  *0^  (6300A),  *0g 
(5577A)  and  vibrational  states  of  N2)  the  escape  rate  of  soft 
electrons,  for  instemce,  to  the  dark  side  of  the  magnetically 
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Figure  10.  Coapariaon  of  Energy-Treuiafer  Rates  for 
Photoelectrona  Interacting  with  Thermal  Electrons. 
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conjugate  ionosphere  is  determined  by  the  ambient  plasma.  In 
fact,  for  high  invariant  latitudes,  the  escaping  electrons 
from  the  lower  ionosphere  (~250  km)  travel  through  either  the 
trough  or  the  plasmasphere,  where  the  main  interaction  is 
electron-plasma,  before  reaching  the  conjugate  dark  side.  The 
part  of  the  study  here  was  to  develop  a realistic  model  of 
such  an  interaction  and  incorporate  that  in  the  transport  cal- 
culations. To  this  end,  it  was  discovered  by  this  author  that 
the  electron  transport  equation  (see  Eq.  (17))  could  be  solved 
(and  was  solved)  by  a variational  principle  combined  with  the 
Green's  function  approach.  The  angular  differential  cross- 
sections  could  then -be  incorporated  as  an  empirical  fit  to  the 
data  points.  The  use  of  the  Green's  function  was  made  to  set 
up  the  boundary-value  problems  for  the  three  regions  of  the 
upper  atmosphere.  Another  point  of  importance  which  was 
deemed  absolutely  necessary,  was,  to  account  for  the  transport 
of  electrons  in  the  regions  of  large  magnetic  dip  angles  (for 
instance  the  auroral  oval  and  the  polar  region).  This  was  due 
to  the  fact  that  the  atmospheric  density  profile  changed  as  a 
function  of  altitude  and  thus  presented  an  inhomogeneous  in- 
teracting medium.  Thus  far,  we  have  been  able  to  account  for 
the  inhomogeneity  only  under  the  approximation  that  within 
about  two  mean  free  paths  of  electrons  the  fractional  abun- 
dance of  various  species  can  be  approximated  by  a local  mean 
value.  For  altitudes  above  the  perigee  of  INJUNS  satellite 
(-660  km)  this  approximation  is  good.  Between  250-600  km  is 
the  transition  region  where  a great  deal  of  care  is  required 
to  accurately  calculate  the  collision  rates  of  hyperthermal 
electrons  as  a function  of  their  releucation  lengths.  Below 
250  km,  the  so  called  "local  hypothesis"  is  fairly  valid, 
i.e.  the  hyperthermal  electzK>ns  are  essentially  thermalized 
in  the  localized  region  so  that  the  transport  effects  can  be 
neglected.  As  a final  point  regarding  our  task,  electric  and 
magnetic  fields  were  incorporated  in  the  calculations  (see 
Eq.  (19)).  By  the  very  nature  of  the  Lorentz  force  term,  the 
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resulting  £ x B drifts  and  the  loss  of  the  azimuthal  symmetry 
of  the  spiralling  motion  of  electrons  makes  the  problem  of 
multi-dimensions  both  in  the  spatial  euid  the  momentum  coordi- 
nates. We  have  succeeded  to  some  extent  in  simplifying  the 
transport  equation  by  various  transformations.  However,  a 
complete  theory  is  not  available  yet. 

To  summarize  our  contributions  to  the  AF6L  project 
on  the  morphology  of  the  electron  transport,  our  performance 
for  the  net  period  of  foxir  months  is  as  follows: 

1)  For  a linearly  anisotropic  pitch  angle  scattering 
we  have  a complete  set  of  transport  models  for  the 
temperate,  the  precipitation  and  the  polar  region. 
Of  these  the  computer  program  for  the  temperate 
zone  has  been  given  to  the  Electrical  Processes 
Branch  (R.  C.  Sagalyn)  of  AFGL.  The  program  needs 
to  be  debugged.  However,  two  subroutines  in  that 
progreun,  which  calculate  the  electron  impact  sec- 
ondary ionization  and  the  continuous  rate  of  energy 
loss  of  hyperthermal  electrons  to  the  cool  ambient 
plasma,  are  debugged  and  working. 

2)  A discovery  was  made  that  a variational  principle 
can  be  used  (and  was  used)  to  treat  the  realistic 
form  of  the  pitch  angle  scattering  in  the  transport 
calculations.  A paper  is  being  written  on  that  ap- 
plication and  will  be  submitted  to  a suitable 
journal. 

3)  Study  has  been  initiated  to  incorporate  electric 
and  magnetic  fields  in  the  transport  phenomenon. 

4)  In  cooperation  with  R.  C.  Sagalyn  of  AFGL,  Drs. 
William  Burke  and  Lalita  Rao  (of  Regis  College) , 
INJUNS  data  was  analyzed  for  the  study  of  SAR  arcs 
and  the  electron  heating  in  the  conjugate  dark  side 
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of  the  ionosphere  by  the  photoelectrons  transported 
from  the  sunlit  side.  Abstracts  of  two  papers  sub- 
mitted to  the  J.  Geophys.  Res.  (1977)  are  given 
below. 
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INJUN  5 LOW>£NERGY  PLASMA  OBSERVATIONS 
DURING  A MAJOR  MAGNETIC  STORM 


Electron  densities  and  temperatures  as  well  as  the 
omnidirectional  flux  of  positive  ions  with  E > 28  eV  were  mea- 
sured by  the  spherical  Langmuir  probes  aboard  Injun  5 at  alti- 
tudes greater  than  2000  km  during  the  October-November  1968 
geomagnetic  storm  period.  During  the  early  phases  of  the 
storm,  the  electron  density  in  the  trough  decreased  and  the 
temperature  increased.  As  the  storm  progressed,  the  position 
of  the  trough  nioved  equatorward.  Plasma  erosion  was  observed 
to  the  invariant  latitude  of  40®  during  the  early  recovery 
phase.  The  latitude  of  the  transition  between  light  and  heavy 
ion  dominance  also  moved  equatorward,  but  recovered  at  a 
slower  rate  than  the  position  of  the  electron  trough.  Most  of 
the  hyperthermal  ions  measured  near  the  trough  were  due  to 
ring-current  particles  reaching  to  the  satellite's  altitude. 

The  minimum  electron  densities  in  the  trough  were  measured  to 
be  within  1®  of  latitude  of  the  maximum  ion  flux.  The  maximum 
electron  temperatures  were  observed  several  degrees  equatorward 
of  the  maximum  ion  flux.  At  the  reported  time  and  latitude 
of  an  SAR  arc,  an  electron  temperature  of  about  4700®K  was 
observed,  whereas  in  the  absence  of  the  ring-current,  a tem- 
perature of  about  2000®K  would  be  expected.  The  observations 
^uc'e  also  used  to  evaluate  a method  for  calculating  the  posi- 
tion of  ring  current  using  magnetic  fluctuations  observed  at 
ground  level. 
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THERMAL  AND  HYPERTHERMAL  ELECTRON  DISTRIBUTIONS  IN  THE 


MIDNIGHT  SECTOR  OF  THE  WINTER  TOPSIDE  IONOSPHERE 


A gridded  spherical  electrostatic  analyzer  aboard 
Injun  5 has  been  used  to  measure  fluxes  of  thermal  and  hyper- 
thermal  electrons  at  sub-auroral  altitudes  in  the  midnight 
sector  of  the  northern  ionosphere  between  altitudes  of  2500 
and  850  km.  Hyperthermal  fluxes,  consisting  of  energetic 
photo electro ns  that  have  escaped  from  the  sunlit  southern 
hemisphere  are  observed  along  orbits  over  the  Atlantic  Ocean 
and  North  America  but  not  over  Asia.  The  electron  tempera- 
tures near  2500  km  have  their  highest  values  at  trough  for 
all  longitudes.  In  the  longitude  sector  to  which  conjugate 
photo electro ns  have  access,  T^  - 4000®K  at  2500  km  and  ~3000®K 
at  1000  km.  For  regions  with  the  conjugate  point  in  darkness 
T^  = 2300®K  over  the  1000-2500  km  altitude  range.  Effective 
spectral  characteristics  of  the  photoelectrons  are  studied  as 
functions  of  latitude  and  altitude.  Based  on  these  observa- 
tions, it  is  concluded  that:  Cl)  Conjugate  photoelecti^ons 
are  not  the  major  contributors  to  trough  heating;  (2)  Heat 
conduction  rather  than  local  heating  by  conjugate  photoelec- 
trons is  responsible  for  electron  temperature  distributions 
observed  in  regions  with  sunlit  conjugate  points. 


T 
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